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Abstract. Wc update constraints on cosmic opacity by combining recent SN Type 
la data compilation with the latest measurements of the Hubble expansion at redshifts 
between and 2. The new constraint on the parameter e parametrising deviations 
from the luminosity-angular diameter distance relation (d^ = eJa(1 + z) 2+e ), is 
e = —0.04+0 07 (2-f)- For the rcdshift range between 0.2 and 0.35 this corresponds to 
an opacity At < 0.012 (95% C.L.), a factor of 2 stronger than the previous constraint. 
Various models of beyond the standard model physics that predict violation of photon 
number conservation contribute to the opacity and can be equally constrained. In 
this paper we put new limits on axion-like particles, including chameleons, and mini- 
charged particles. 



1. Introduction 

Cosmological observations provide constraints on different distance measures: 
luminosity distance (as provided e.g., by supernovae), angular diameter distance (as 
provided e.g., by baryon acoustic oscillations) and even on the expansion rate or 
the Hubble parameter as a function of redshift z. Both luminosity distance and 
angular diameter distance are functions of the Hubble parameter. While combining 
these measurements helps to break parameter degeneracies and constrain cosmological 
parameters, comparing them helps to constrain possible deviations from the assumptions 
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underlying the standard cosmological model (e.g. isotropy), or to directly constrain 
physics beyond the standard model of particle physics (e.g. couplings of photons to 
scalar or pseudo-scalar matter). 

The Etherington relation [1] implies that, in a cosmology based on a metric theory 
of gravity, distance measures are unique: the luminosity distance is (1 + z) 2 times the 
angular diameter distance. This is valid in any cosmological background where photons 
travel on null geodesies and where, crucially, photon number is conserved. 

There are several scenarios in which the Etherington relation would be violated: for 
instance we can have deviations from a metric theory of gravity, photons not traveling 
along unique null geodesies, variations of fundamental constants, etc. In this paper we 
want to restrict our attention on violations of the Etherington relation arising from the 
violation of photon conservation. 

A change in the photon flux during propagation towards the Earth will affect 
the Supernovae (SNe) luminosity distance measures but not the determinations of the 
angular diameter distance. Photon conservation can be violated by simple astrophysical 
effects or by exotic physics. Amongst the former we find, for instance, attenuation 
due to interstellar dust, gas and/or plasmas. Most known sources of attenuation are 
expected to be clustered and can be typically constrained down to the 0.1% level [2j[3]. 
Unclustered sources of attenuation are however much more difficult to constrain. For 
example, gray dust [I] has been invoked to explain the observed dimming of Type la 
Supernovae without resorting to cosmic acceleration. 

More exotic sources of photon conservation violation involve a coupling of photons 
to particles beyond the standard model of particle physics. Such couplings would 
mean that, while passing through the intergalactic medium, a photon could disappear 
-or even (re) appear!- interacting with such exotic particles, modifying the apparent 
luminosity of sources. Here we consider the mixing of photons with scalars, known as 
axion-like particles, and the possibility of mini-charged particles which have a tiny, and 
unquantised electric charge. A recent review [5] highlights the rich phenomenology of 
these weakly- interacting-sub-eV-particles (WISPs), whose effects have been searched for 
in a number of laboratory experiments and astronomical observations. In particular, 
the implications of this particles on the SN luminosity have been described in a number 
of publications [71 El El HQ] . 

One of the most interesting features of these models is that the exotic opacity 
involved could in principle "mimic" the value of a non-zero cosmological constant 
inferred from SNe measurements. However, this possibility can already be excluded 
(at least in the simplest WISP models) by the absence of distortions in the CMB or the 
spectra of quasars for axion-like-particles, and by arguments of stellar evolution in the 
case of mini-charged particles. 

In this paper we use improved bounds on cosmic opacity to further constrain the 
existence of exotic particles which can couple to the photon. The rest of the paper 
is organised as follows. In section 2 we update constraints on transparency from the 
latest available data. In section 3 we discuss the implications of this for axion-like 
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particles and chameleons, and in section 4 we consider mini- charged particles. We then 
forecast, in section 5, how the constraints will improve with distance measures from 
future, planned and proposed, surveys. We conclude in section 6. Sections 3 and 4 
discuss in detail the motivation, modelling and regime of applicability of the beyond the 
standard model physics we consider. Readers with a more focused interest on cosmology 
may concentrate on the beginning of section 3, sub-sections 3.4, 3.5 and figures 2, 3, 4, 
5, 6. Appendix A summarises the cosmologically-relevant results of sections 3 and 4. 
Precursors of this paper can be found in [TTJ [12], [13], [HJ EES] . 

2. An update on cosmic opacity constraints 

In reference [15], the authors use Type la SN brightness data (namely the SCP Union 
2008 compilation [16]) in combination with measurements of cosmic expansion H(z) 
from differential aging of luminous red galaxies (LRGs) [TTJ [18] to obtain constraints on 
non-trivial opacity, at cosmological scales. The basic idea is to study possible violations 
from the "Etherington relation" [T], the distance duality between luminosity distance, 
<ii, and angular diameter distance, <i^: 

d L (z) = (1 + zfd A {z) . (1) 

This identity depends only on photon number conservation and local Lorentz invariance. 
It holds for general metric theories of gravity, where photons travel along unique 
null geodesies. Since Lorentz violation is strongly constrained for the low energies 
corresponding to optical observations [19], the study of possible violations of Eq. (TTJ) 
through SN observations directly constrains photon number violation. Any such 
systematic violations can then be interpreted as an opacity effect in the observed 
luminosity distance, parametrised through a generic opacity parameter, t(z), as: 

^ 2 L,obs = d 2 Ljrue eT ^ • (2) 

Note that our "opacity" can have in principle both signs. In other words, this 
parametrisation also allows for apparent brightening of light sources, as would be the 
case, for example, if exotic particles were also emitted from the source and converted 
into photons along the line of sight [9]. From Eq. (j2j) it is clear that the inferred distance 
moduli for the observed SNe picks an extra term which is linear in r(z): 

DM obs {z) = DM true (z) + 2.5[loge]r(z) . (3) 

On the other hand, one can also use other determinations of distance measures, 
which are independent of r, to constrain possible deviations from Eq. ([[]). This approach 
was initiated in reference [H] (see also [TTJ [121 I20J [21] for related earlier work) where 
the authors used measurements [22] of the baryon acoustic oscillation (BAO) scale at 
two redshifts, namely z = 0.20 and z = 0.35, to obtain a parameterization-independent 
upper-bound for the difference in opacity between these two redshifts, At < 0.13 at 
95% confidence. In reference [15] this constraint was improved (and also extended over 
a wider redshift range, but for a general parameterised form for r) by using, instead 
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of measurements of the BAO scale at these two redshifts, measurements of cosmic 
expansion H(z) from differential aging of LRGs at redshifts z < 2. This method of 
distance determination relies on the detailed shapes of galaxy spectra but not on galaxy 
luminosities, so it is independent of r. 

In particular, the authors introduced a parameter e to study deviations from the 
Etherington relation of the form: 

d L {z) = d A (z)(l + z) 2+ % (4) 

and constrained this parameter to be e = — 0.0li^09 a ^ 95% confidence. Restricted to the 
redshift range 0.2 < z < 0.35, where t(z) = 2ez + 0(ez 2 ), this corresponds to At < 0.02 
at 95% confidence. In the following sections, we will apply similar constraints on different 
parametrisations of r which correspond to particular models of exotic matter-photon 
coupling, namely axion-like particles (ALPs), chameleons, and mini-charged particles 
(MCPs). 

Before moving to these models, we briefly update the above constraint on e using 
the latest H(z) data [23], which include two extra data points at redshifts z = 0.48 
and z = 0.9, as well as the latest determination of H Q [23]. Even though the addition 
of these two extra data points alone significantly improves the constraints of reference 
|15] . the effect of Ho is also quite significant, because it acts as an overall scale in the 
distance measures, which is marginalised over a Gaussian prior, and the measurement 
error in this determination is about half of that of the HST Key Project determination 
[25] used in [T5j . 

Fig. [T] shows the updated constraints obtained using the above data in combination 
with the SCP Union 2008 Compilation [16] of type la Supernova dataljl, compared to the 
previous constraints of reference [15]. On the left, the darker blue contours correspond to 
the (two-parameter) 68% and 95% joint confidence levels obtained from SN data alone, 
while lighter blue contours are the corresponding confidence levels for H(z) data. Solid- 
line transparent contours are for joint S~N+H(z) data. For comparison we also show the 
previous H(z) and SN+if(,?:) contours in dotted and dashed lines respectively. On the 
right we show one-parameter (marginalized over all other parameters) constraints on e, 
again for the current analysis (solid line) and for that of reference [15] (dotted). For the 
reader familiar with Bayesian methods, this plot corresponds to the posterior 

P(e|S,E)= f I P(n m ,H \E)P(e,n m ,H \S)dn m dH , (5) 

J Q m J Ho 

where P(Q m ,H \E) and P(e, Q m , H \S) are the posterior probabilities for the 
corresponding model parameters, given the H(z) (Expansion) and SN (Supernovae) 
data respectively. These are given by the likelihoods of the two data sets in the model 
parameters, assuming Gaussian errors and using flat priors on all three parameters. In 
particular, we have taken e E [-0.5, 0.5], Vt rn E [0, 1] and H E [74.2-3x3.6, 74.2+3x3.6] 

X Note that we decide to use the SCP Union 2008 supernova sample rather than the more recent 
SDSSII sample [26], because the Union sample extends to higher redshift and is thus best suited to be 
combined with the H(z) data. 
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(Riess et. al. [21]), all spaced uniformly over the relevant intervals, in a flat ACDM 
model. Similarly, the solid line transparent contours on the left plot of Fig. [1] correspond 
to taking only the integral over Hq in the right hand side of Eq. (jHJ), yielding, therefore, 
the two-parameter posterior P(e, f2 m |S, E). 
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Figure 1. Updated constraints of reference [15], using the latest H{z) data [23] and 
the Riess et al. determination of Hq [24] in combination with the SCP Union 2008 SN 
la compilation. Left: Two-parameter constraints on the e — Q m plane. Darker blue 
contours correspond to 68% and 95% confidence levels obtained from SN data alone, 
lighter blue contours are for H(z) data, and solid line transparent contours are for 
joint SN+if (z). Previous H(z) and joint SN+i?(z) from [T5] are shown in dashed and 
dotted lines respectively. Right: One-parameter joint constraints on e for the current 
analysis (solid line) and that of reference [T5] (dotted line) . The dashed line shows the 
95% confidence level, A^ 2 = 2. 



As seen in Fig. [IJ the improvement in these constraints is significant. The new 
result on e, marginalised over all other parameters, is e = — 0.04+q;q7 at 95% confidence, 
which for redshifts between 0.2 and 0.35 (currently probed by BAO data), corresponds 
to a transparency (i.e., r > 0) bound At < 0.012, a factor of two tighter than the result 
in reference |15)§l . We now move on to study more general parametrisations of cosmic 
opacity, tailored for particular models of exotic matter coupled to photons. 

3. Axion-like Particles and Chameleons 

New scalar or pseudo scalar particles from physics beyond the standard model, here 
denoted as 0, may couple to photons through 

Cscalar = 7^ V> (6) 

§ Note that the data slightly favour negative e (thus the much stronger constraint on a positive At), 
but only at < l-cr level. 
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where M is the energy scale of the coupling (another widely used notation is = 
1/M), F^y the electromagnetic field strength and e^xp the Levi-Civita symbol in four 
dimensions. Such fields are collectively known as Axion-Like Particles (ALPs), as a 
coupling of the form (JTj) arises for the axion introduced by Peccei and Quinn (PQ) to 
solve the strong CP problem [27] . Interestingly, these fields also arise naturally in string 
theory (for a review see [28]). 

Axions, or axion-like-particles, can arise from field theoretic extensions of the 
standard model as Goldstone bosons when a global shift symmetry, present in the high 
energy sector, is spontaneously broken. In the PQ axion case, this symmetry is colour 
anomalous and the explicit breaking makes the axion pick up a small mass. This mass 
is, up to a model-independent constant, proportional to the coupling (JTJ). For a generic 
ALP, however, the mass is in principle independent of the strength of its coupling, and 
in particular can be zero if the releted shift symmetry remains intact. That is, for 
instance, the case of Arions [29], the orthogonal combination of the PQ axion, if there 
are two independent colour anomalous shift symmetries. 

Chameleon scalar fields are another very interesting type of ALPs [30]. They were 
originally invoked in [3TJ [32] to explain the current accelerated expansion of the Universe 
with a quintessence field which can couple to matter without giving rise to large fifth 
forces or unacceptable violations of the weak equivalence principle. The chameleon 
achieves this because its mass depends on the local energy density. The environmental 
dependence of the mass of the chameleon means that it avoids many of the constraints 
on the strength of the coupling, which normally apply to standard scalar and pseudo- 
scalar fields as they are derived from physics in dense environments. For a more detailed 
discussion see [33]. The cosomology of the chameleon was explored in detail in [M], the 
possibility of the chameleon coupling to photons was first discussed in [35] and such a 
coupling was shown to be generic in [30] . 

The Lagrangian terms given above mean that ALPs can affect the propagation of 
photons; in particular, if photons traverse a magnetic field there is a non-zero probability 
that they will oscillate into ALPs [36]. Notice however that only photons polarized 
perpendicular (parallel) to the magnetic field mix with scalar (pseudo-scalar) particles. 
Therefore, the interactions between photons and ALPs in the presence of a magnetic field 
not only imply that photon number is not conserved, but can also alter the polarization 
of the light beam. Both effects have been exploited in many searches for ALPs both in 
the laboratory and in astronomical observations, see [5] for a recent review. 

3.1. Modelling the effects of ALPs 

The presence of ALPs will have an impact on observations of SNe if they are observed 
through intergalactic magnetic fields. In particular, it will lead to changes in the 
observed SN luminosities, in a redshift-dependent way. Many different mechanisms have 
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been proposed which give rise to intergalactic magnetic fields, however we do not yet have 
convincing evidence from observations that they exist. A magnetic field coherent over 
the whole Hubble volume is limited, by observations of the CMB and Faraday rotation, 
to B < 10~ 9 G [37J [38] - Fields with shorter coherence lengths are also constrained. 
In particular, fields coherent on scales ~ 50 Mpc must satisfy B < 6 x 10~ 9 G, while 
fields coherent on scales ~ Mpc must satisfy B < 10~ 8 G [38] • To explain the origin of 
galactic magnetic fields it is expected that intergalactic magnetic fields with coherence 
lengths ~ Mpc are needed [39|. 

In a constant, coherent magnetic field the probability of a suitably polarized photon 
converting into an ALP after traveling a distance L is given by [36J: 

P 7 ^(L) = sin 2 (2^)sin 2 (^) , (8) 

where 

A(L) = (9) 

2Bui 

tan 20 = -—^, (10) 

where u = Itxv is the photon energy, B is the strength of the magnetic field and 
m cff = l m | — k'pl: w hh the mass of the ALP and up = 47r 2 an e /m e the plasma 
frequency of the medium which acts as an effective mass for the photons (a is the fine 
structure constant, n e the local number density of electrons, and m e the mass of the 
electron). In what follows we restrict our attention to very light fields, < up, where 
observations of the opacity of the universe have the most power to constrain the strength 
of the coupling for the ALP to photons. 

However, the intergalactic magnetic field is not coherent from the earth to the 
supernovae. We model its fluctuations using the cell magnetic field model, whereby we 
assume that the magnetic field is split up into a large number of equally sized domains. 
Within each domain the magnetic field is constant, its strength being the same in all 
domains, while its orientation is chosen randomly. The cell magnetic field model is the 
simplest choice we can make to approximate the structure of astrophysical magnetic 
fields, and is commonly used both in the study of ALPs and of astrophysical magnetic 
fields. A more accurate choice for a model of the magnetic field would be to assume a 
power spectrum of its fluctuations. However, at high frequencies, A <C ir/2, the cell and 
power spectrum models give the same results, and at lower frequencies the cell model 
captures all the qualitative features of ALP-photon interactions, but underestimates the 
probability of conversion jlOJ HI]. Therefore, using the cell magnetic field model will 
give rise to conservative constraints. 

Clearly, we also need to know the plasma frequency of the intergalactic medium 
along the line of sight. This is quite a complicated issue because no measurements 
of the electron density are available in the large voids of the interstellar medium. A 
large-scale average can be easily inferred from the total amount of electrons determined 
by CMB estimation of the baryon to photon ratio giving up ~ 1.8 x 10~ 14 eV today, 
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see e.g. [32]. Note, however, that average values up to a factor of 15 smaller were 
considered plausible in [32]. Since there is no easy way out of this conundrum, the 
accepted approach is to assume that u 2 P is homogeneous and equal to the average value. 
To check the dependence of the results on this assumption, we will finally allow a range 
of a couple of orders of magnitude around the average. 

As we are interested in the transparency of the Universe out to redshifts z ~ 0(1) 
we must also take into account the redshift evolution of the environment that causes 
mixing between photons and scalars. Assuming the magnetic fields are frozen into 
the plasma, their strength scales as B(z) = B (l + z) 2 [39] while cu(z) = cu (l + z) and 
ojp(z) = ojp (l + z) 3 (since it is proportional to the electron density). Here, the subscript 
indicates values in the present epoch. The physical length of a magnetic domain scales 
as L(z) = Lq(1 + z) _1 as long as it is smaller than the Hubble radius. Then the two 
parameters that appear in the probability of conversion (jSJ), 9 and A, redshift as 

A(*) =A (l + z), (11) 
tan2#(^) = tan 2^ . (12) 

There are two limits in which the expression for the conversion probability in 
one domain simplifies notably, the coherent and incoherent regimes. In the coherent 
regime, the argument of the sinus in Eq.® is smaller than 1 and, taking sinx ~ x, the 
probability takes a very simple expression: 

W«) « (^jr) 2 (! + z ? (coherent) . (13) 

On the other hand, if the argument of the sinus - which is energy dependent - is very 
large, then a large number of oscillations will happen in a finite energy bin, which would 
average out the sinus to 1/2. In this case we find 



2 



P-^Az) « - „ (incoherent) . (14) 



2 



\Mu 2 Pfl ) 



These approximations are only valid for small values of P^^. Finally, note that from 
now on we will drop the subscript for today's values of the various parameters and 
make the redshift dependence explicit. 

In the above limits, the redshift dependence is very simple and a system of 
axion-like particles and photons can be evolved analytically through a large number 
of randomly oriented magnetic domains. Let us introduce the notation P for the 
transition probability P 7 _>0 in one domain today. Assuming that the magnetic fields and 
intergalactic medium do not evolve with redshift, the photon survival probability was 
first computed in [33] • This is then exactly valid in the incoherent regime of Eq. ( TT4l) . 
The finite probability of conversion gives rise to an apparent change in luminosity 
distance. In particular, if photons are converted to ALPs along the line of sight, then 
the inferred and true and luminosity distance squared (cf Eq. (J2J)) will differ by a factor 
V(z), which in this case reads: 

V(z) = A + (1 - A) exp , (15) 
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where y(z) is the comoving distance to the source. Physically, Viz) is the average 
probability that a photon emitted by a supernovae at redshift z is observed by us after 
traversing the magnetic fields in the intergalactic medium. The above formula is valid for 
small P; in the case where P is of order unity, one should replace 3P/2 — > — ln(l— 3P/2). 
We have allowed for an initial flux of axions I<j>(zi) and definecljj] 



■■ A -A i+ m)- (i7) 

In the coherent regime, we wish to allow for additional effects due to the evolution 
of the magnetic fields with redshift. The fluxes of scalars and photons at the end of the 
n-th domain are related to the fluxes at the beginning of the domain by 

I 1 {z)\ (l-\P(l + z? P(l + zf \(l 1 {z)\ 

I ^< 1+2 > 2 '-'I'+'mwj,,, 

and we want to multiply a large number of these matrices together. Usefully, the matrix 
is diagonalisable by a redshift independent transformation, so that, after passing through 
./V domains, the photon survival probability can be easily shown to be 

N / a -\ 

(19) 



V(z) = A+(l-A)f[(l-~P(l + z j )^ 



where {zj} is a collection of redshifts in the range (0,2:), equally spaced in comoving 
length. If the number of domains is large, we can approximate the product by an integral 
to get 

V{z) = A + (1 - A) exp (-jfj H ^~ H0 ) ■ (20) 

This resembles the expression for the photon survival probability when the evolution of 
the background is neglected, equation (fl5l) . but has a stronger ^-dependence at large 

| Very recently it was pointed out that this formula is actually an averaged formula over different 
realisations of the configuration of magnetic field domains along a line of sight [15J 125] . As such, it 
is in principle not valid for a single source, whose light only travels through a concrete realisation 
of the magnetic field domain structure. The authors of [J5] calculated an analytical estimate for the 
dispersion around the mean - although the distribution is non-Gaussian and in general asymmetrical 
so a meaningful confidence interval has to be computed numerically - for A = 2/3. Using x — Py(z)/L, 
the result reads 

SP = y[R=T» with ^49 + 50cxpM,/2) + 6cxpH^ ^ 

luo 

which falls in the range ~ 

(0,l/(3\/B)), i-e. the dispersion is below 15%. In the Hubble diagram, this 
corresponds to a maximum dispersion of 0.2 magnitudes, which is the typical dispersion in the SN 
data wc use. This otherwise suggestive fact, implies that confronting the observed with the predicted 
dispersion can only potentially constrain a regime of large values of x - what we will later introduce as 
strong mixing regime - and most likely with small significance. We conclude that it is more promising 
to confront the mean value, given by Eq. (|15[) . with observations. In fact, as we shall see below, when 
constraining "opacity" as a function of redshift, for every reasonable redshift interval we may consider, 
there will be several data points corresponding to different supernovae at different positions in the 
sky, yielding an effective average; therefore in the following we can obviate the implications of the 
uncertainty in the value of 8P. 
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redshifts. Both formulae give the same results in the small z regime, since y(z) ~ z/Hq 
and (H/H - l)/fi m ~ 3z/2. 

The redshift evolution of Eqs. ( 1151) and ( 1201) is absent when A = 1. When this is 
the case, the initial flux of photons and scalars is already thermalised I^{zi) = 2I^(zi) 
so that, on average, the effect of photons converting into scalars is compensated by that 
of scalars converting into photons, and thus no net effect of the mixing is seen. 

In summary, the luminosity distance to a supernova is modified by an overall, 
redshift-dependent factor: 

d L ^d L /^Vjzj. (21) 

As mentioned in section [21 the effect of ALP-photon mixing described above can be 
interpreted as "opacity" , generally of both signs, so that 

t(z) = -hxV(z) . (22) 

We now proceed to constrain the possibility of mixing between scalars and photons, 
through their effects on cosmic opacity. Again, there are two regimes in which analytical 
insight can be reasonably expected. 



3.2. Weak mixing 

We begin by considering the case where the sum of the 7—7-0 conversion probabilities 
in all the domains is smaller than unity (iVP 7 ^ < 1). In this limit the effects of 
ALP-photon mixing are always small so it is known as weak mixing regime. When the 
redshift evolution of the background is neglected, the probability of photon survival can 
be found [331 SOI E] to be 

P(z) = l-(l-A)l(y&)P. (23) 



2 V L 

In the low redshift regime, both the redshift dependent and independent equations give 
the same opacity 

t{z) = (1 ~ a) \4^l z ■ (24) 

Note that this is exactly of the form used in section 2, t(z) = 2ez. We can already 
obtain a first estimate of our bounds by using our improved constraint e < 0.04 at 95% 
confidence. Using H ~ 74.2 km/s Mpc -1 we obtain 

P < 4 x 10- 5 -^ (z^t) ■ (25) 
Mpc \1- AJ ( J 

Note that A ~ 1 cannot be constrained. As mentioned before, this situation 
corresponds to the initial flux of axions and photons having almost thermalized 
abundances (see Eq. (I17jl ). A thermalized, or nearly thermalized, axion/photon flux 
has very small redshift dependence (since the mixing tends to thermalization, but this 
is almost complete before leaving the SN source) and thus distance measures have no 
constraining power. 
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Two further comments are in order. Our bound of Eq. (1251) corresponds to the 
argument of the exponential taking values around — 0.081z/(l — A). This corresponds 
to the start of the exponential regime for the larger redshifts, so it is consistent with the 
Taylor expansion adopted, unless A is fine-tuned to 1. Finally, note that the bound on 
e has been obtained considering redshifts up to z ~ 1.5, while the validity of Eq. (121]) is 
only ensured for small redshifts. For this reason this bound should be considered as an 
order-of-magnitude estimation; we present an accurate, numerical study, in section 13.41 



3.3. Strong Mixing 

The other analytically analysable limit of Eqs. ( jT5"l) and ( 120]) is when the sum of the 
conversion probabilities in all domains is very large, NP 3> 1. If the mixing between 
axion-like particles and photons is strong, then, on average, one-third of any initial flux 
will be converted into axion-like particles and two thirds into photons. This can cause 
substantial changes to the apparent supernovae luminosities. It was initially speculated 
[6] that this could account for the dimming of SNe without the existence of dark energy, 
however this has now been excluded by observations of the CMB [UJ . If there is a large 
initial flux of axion-like particles from SNe, then the SN images can be brightened by 
strong mixing between ALPs and photons in the intergalactic medium [9]. 

In the strong mixing limit both equations ( fl~5]) and ( 1201 give the same result, V ~ A, 
that is, a constant independent of redshift. It is clear that our methods cannot constrain 
such a possibility, which would be equivalent to a change of the normalization of the SN 
luminosity. Imposing NP > 1 at the closest SN redshifts up to z ~ 0.01, we find 

P> 0.015-^-. (26) 
Mpc v 1 

We cannot constrain this range of probabilities as the mixing is too strong and the 
redshift dependence of the opacity is washed out. 



3.4- Numerical results 

In this subsection we turn our estimates into real constraints by means of a full likelihood 
analysis. Note that our results only depend on the conversion probability per comoving 
length P/L, so we cannot constrain P and L independently. 

In Fig. [2] we show our results marginalized over Q m and Hq. The dark and light 
contours represent 68% and 95% joint confidence levels respectively, using the SN data 
only (left) and joint SN+if (z) data (right). In the upper panels, we have used Eq. ( 1201) . 
thus taking into account the redshift dependence of the background, while in the lower 
panels we used Eq. ( fl~5]) . ignoring these effects. 

The strong and weak mixing limits described in sections 13.31 and 13.21 are clearly 
visible in Fig. [21 We represent them schematically by the green and red lines below and 
above the boundaries of our plots respectively. For very small conversion probability 
P < 10~ 5 we are in the weak mixing limit and our constraints become weak because 
of the lack of photon-ALP mixing. On the other extreme, when the probability is 
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Figure 2. Two parameter constraints in the A — P/L plane for general ALPs 
considering rcdshift dependence of the background (upper panels) and neglecting it 
(lower panels). Contours represent the 68% (dark) and 95% (light) confidence levels. 
For the left panels we have used SN data only, while in the right panels we show 
the joint SN+H(z) analysis. The contours are marginalized over cosmologies and Hq. 
While rcdshift dependence introduces 0(1) effects, which change the structure of the 
68% CL contours, the 95% CL contours are very similar in the two cases. 

very strong (P > few x 10 -2 ) the photons and axions mix until thermalization and 
the redshift dependence of the opacity is lost, so, again, our constraints become weak. 
Finally we can also observe that a band around A = 1 (yellow line) again cannot be 
constrained. A = 1 means an equilibrated photon-ALP flux from the SNe such that 
photon— ALP conversions are compensated with the reverse process, making the photon 
number constant, i.e. redshift-independent. 

The first notable feature is that, as expected from our estimations in the first 
section, our 2-a bounds (white regions) are significantly improved when including 
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the H(z) data (right panels) compared to using SN data alone (left panels). The 
improvement is particularly visible in the weak mixing regime. This will show more 
clearly in the constraints on the P/L — Q m plane, which we will present below. We shall 
comment on this improvement separately in each case of interest. 

The second effect is the importance of including the redshift dependence of the 
magnetic fields, the photon frequency and the physical length of domains, that is, the 
difference between using Eqs. ffl5|) and ( 120]) . The 95% C.L. constraints are improved only 
slightly by including the redshift dependence (upper panels), compared to a background 
that does not evolve (lower panels), the exception being the factor of 3 improvement at 
A ~ 2/3 in the SN+if (z) analysis (right panels). 

A peculiar feature appears in the upper left plot, corresponding to a redshift- 
dependent background and using SN data only. There, the 68% C.L. contour shows 
not one, but three regions where photon-ALP mixing improves the fitting of the 
data compared to a standard ACDM cosmology. The lower-left region corresponds 
to the parameters invoked by Csaki et al. [6] to explain SN dimming without cosmic 
acceleration. It is remarkable that the data show this preference much more sharply in 
the redshift-dependent background case given by Eq. |20] than in the simpler redshift- 
independent case used in [6] (Eq. [T5|) . As we show in the next section, our numerical 
analysis implies that this island in parameter space corresponds to Q m < 1, which, in 
the absence of a cosmological constant, implies a non-flat geometry. Most importantly, 
the joint SN+if (z) analysis (upper right panel) rules out this region. 

The preferred region at large A and small P/L also deserves some comments. A 
value of A greater than unity produces an increase of the SN luminosity with redshift 
because the SNe would shine more ALPs than photons. In Fig. [T]we have shown that, 
because of the pronounced degeneracy, SN data alone show a slight preference for this 
scenario since the 68% C.L. contour is almost completely located at negative values of 
our opacity parameter e. In the joint SN+H(z) analysis, the trend is softened but, still, 
small negative values (slight SN brightening) are slightly preferred, even though this is 
not statistically significant. Therefore, in the remaining panels of Fig. [2] this possibility 
cannot be excluded, but the statistical preference for this region decreases compared to 
the rest of the allowed parameter space. The fact that photon-ALP mixing alleviates 
the tension between SN data and standard rulers was already pointed out in [9]. The 
value of P/L ~ 2 x 10~ 5 Mpc -1 implies that the required effect on the opacity is small. 
In the weak mixing limit, the opacity is proportional to (1 — A)P/L (see Eq. ( EH1) ). 
making A and Pj L correlated. This region is therefore expected to extend further right 
and down. Notably, these small values of P/L are not excluded by other arguments 
(see next section), although it seems difficult to conceive a model where A is sufficiently 
large. 
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Figure 3. Confidence levels (68% and 95%) on the P/L — f2 m plane for the simplest 
axion- like-particle model (A = 2/3). The small and large P regions correspond to the 
weak and strong mixing regimes respectively. Dark blue contours show constraints from 
SN data only, light blue from H(z) data, and solid line contours from joint SN+H(z). 
In the left panel, the redshift dependence of the background is taken into account, 
while in the right panel these affects are ignored. 

3. 5. Constraints 

We can now interpret our results in terms of constraints on the physical parameters of 
our Lagrangian. There are several cases to consider. 

3.5.1. Axion-like-particles: If ALPs have no other interactions than the two-photon 
coupling, then they were shown to contribute very little to the SN luminosity [44], 
corresponding to the case A ~ 2/3. It is evident from Fig. [2] that only a range of 
conversion probabilities around P/L ~ (9(10~ 3 ) can be excluded. Fig. [3] shows 1- and 
2-er, two-parameter likelihood contours on the P/L — Q m plane after marginalization 
over Hq. Note that even the SN constraints alone (dark blue contours) rule out this 
model as an alternative to a cosmological constant, at greater than 3-a significance. 
However, there is still significant degeneracy in the weak mixing regime, and a value for 
Q m greater than 0.8 is still allowed at the 2-a level, this constraint being slightly weaker 
when the background redshift dependence is taken into account (left). Including the 
H(z) data (light blue contours) breaks this degeneracy, yielding strong joint constraints 
in P/L — Q m (solid line contours). 

We can now translate the bounds on Pj L into bounds on the strength of the ALP 
coupling to photons. Since the coupling always appears multiplied by the magnetic field 
(which is also unknown) we find convenient to quote bounds on the combination B/M. 
Let us also define appropriately normalised values of the magnetic field strength B and 
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the energy-scale of the the axion-photon coupling M, as 
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In Fig. H] (left) we show our constraints for the case L = 1 Mpc as a function of the 
uncertain value of the average electron density, or, equivalently, the plasma frequency. 
To get rid of the oscillations of Eq. (jSJ), which not only will be averaged out by energy 
binning but also by small fluctuations in the sizes of the domains and the values of 
the plasma frequencies, we propose the substitution sin 2 x — > (1 — exp(— 2x 2 ))/2, which 
reproduces the coherent and incoherent limits. The exclusion limit is a horizontal band 
which bends upwards around n e ~ 0.2 x 10 -7 cm -3 . The horizontal part corresponds 
to the coherent case, where the n e dependence drops out of P, while the diagonal band 
corresponds to the incoherent case where n e suppresses P. Note that the average electron 
density today is n e ~ 10 -7 cm -3 , near the transition between the two regimes. 

In the same figure, we have also reproduced the constraints of Mirizzi et al [47] from 
distortions of the CMB (region above the dashed line) and those of [IE] (blue region) 
from QSO spectra (see also footnote 2 in [47] and [491 [50]). For n e > 10~ 9 cm -3 our 
bounds are stronger than the CMB ones while still competitive with the QSO bounds. 
Our approach provides a complementary, independent way to obtain these constraints. 
Each of the three approaches reported in the figure (especially the present constraint and 
the QSO one) is affected by different, unrelated, systematics: their agreement promotes 
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Figure 4. Bounds for the product of ALP coupling times magnetic field B n o/MiQ as 
a function of the average electron density or corresponding plasma frequency. The gray 
regions are excluded at 95% C.L. for L = 1 Mpc in the redshift-indcpcndcnt (dark) and 
dependent (light) cases. Also shown are the regions constrained from CMB [47] (above 
the dashed line), which dominate at low n e , and QSO |48j (blue) spectral distortions. 
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one's confidence in these results. 

3.5.2. Chameleons: Unlike the simple ALP case studied above, the chameleon 
Lagrangian contains non-linear self-interactions of the scalar field in order that the mass 
of the scalar may become dependent on the density of its environment. This introduces 
the possibility of having A ^ 2/3. 

There are three possible cases: i) either A ~ 2/3 because few chameleons are 
produced in the SN, ii) A — 1 -C 1 because they interact so strongly in the SN that 
photons and chameleons thermalize their fluxes within the SN, and iii) the intermediate 
case where a significant number of scalar particles are produced in the interior of the 
SN, but, yet, photon-chamelon interactions are not strong enough to thermalise the 
chameleon population with that of the photons before they leave the SN. 

The first case is morphologically equivalent to the previous ALP case so the 
conclusions of the last section hold. In the second case we see from Fig. [2] that we 
cannot constrain any value of the probability of conversion. Our sensitivity is at most 
\A — 1| > 0.1 but only in a very narrow range of P around 10 _3 L/Mpc. Therefore, 
we cannot exclude the possibility that photons and chameleons mix strongly in the 
intergalactic medium, even though, in order to realise this scenario fully, more work 
is required to understand the possible fluxes of chameleons from SNe. For the third 
case we can obtain a constraint on the P/L — Q m plane by marginalizing over A in the 
physically reasonable range [2/3,4/3] (refer to Eq. ( fl71) ). Fig. Oshows our results in this 
case together with the corresponding constraint on the A — Q m plane after marginalising 
over P/L in the range [10~ 5 , 10 -1 ]. 

Note that, in agreement with our analysis of the simplest ALP case with A = 2/3, 
our constraints exclude these models as an alternative to a cosmological constant at 
greater than 3-cr if the Universe is set to be spatially flat, while Q m > 0.8 is excluded 
at 2-cr even for SN data alone, and at much higher statistical significance when H(z) 
data are included. This is true both when the redshift dependence of the background 
is included (upper panels) or neglected (lower panels), the constraint in the latter case 
being somewhat stronger. 

4. Mini-Charged Particles/Hidden Photons 

New particles with a small unquantized charge have been investigated in several 
extensions of the standard model [511 [52]. In particular, they arise naturally in 
extensions of the standard model which contain at least one additional U(l) hidden 
sector gauge group [HU [53] . The gauge boson of this additional U(l) is known as 
a hidden photon, and hidden sector particles, charged under the hidden U(l), get an 
induced electric charge proportional to the small mixing angle between the kinetic terms 
of the two photons. In string theory, such hidden U(l)s and the required kinetic mixing 
are a generic feature [EH [55l [56j [571 EH]- Hidden photons are not necessary however to 
explain mini-charged particles, and explicit brane-world scenarios have been constructed 
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Figure 5. 68% and 95% confidence levels (two parameters) on the P/L — fl m (left) 
and A — Q m (right) planes, for chameleons with a prior A € [2/3, 4/3]. In the upper we 
have taken into account the order 0(1) effects arising from the redshift dependance, 
while in the lowers panels we have ignored these effects. Dark blue contours are for SN 
data only, light blue ones for H(z) data, and black transparent contours are for the 
joint SN+iJ (z) analysis. 

[52] where MCPs arise without the need for hidden photons. 

The existence of low-mass MCPs can have a tremendous impact on photon 
propagation over cosmological distances. Photons from a given source can for instance 
pair produce MCPs with CMB photons 7 + 7cmb ~ > ^ + + leading to a new form 
of opacity. However, this process is generally more noticeable for CMB photons rather 
than those of higher energy, both because the CMB spectrum was measured by the 
FIRAS experiment to be a perfect blackbody with a typical accuracy of 10~ 4 , and, also, 
because the cross-section is inversely proportional to the center-of-mass energy. The 
impact of the existence of MCPs for CMB distortions was studied in [59] , where a limit 
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for the minicharge q e < 4 x 10~ 8 (measured in units of the electron's charge) was derived 
for 4-component Dirac MCPs. 

A more relevant source of opacity was pointed out in [10], following the work of 
[601 [61]. Photons propagating in a background magnetic field can actually pair-produce 
MCPs without the need for a second photon in the initial state. This is due to the 
fact that in a background field energy-momentum conservation is non-trivial. Indeed, 
the magnetic field acts as a refractive medium where both the mini-charged particles 
and photons, acquire a non-trivial dispersion relation, i.e. effective masses. In the most 
interesting case, the effective photon mass is larger than that of a MCP pair, and the 
7 — > ipifj process happens at a rate 

8/3 /2R 2 o 8 \ 1/3 
"= ,wnwm 4±) /, (28) 



«■ (|) r (f ) V 3. 

where q e is the MCP electric charge in units of the electron's charge e, and / is an order 
one factor which depends on the nature of the MCP and the photon polarization with 
respect to the magnetic field, assumed again to be transverse to the photon direction 
of motion (/ = 1,2/3 for parallel and perpendicular polarizations respectively if the 
MCP is a Dirac spinor, and / = 1/12, 1/4 if the MCP has spin-0). T denotes the usual 
T-function. The above formula is valid in the deep non-perturbative regime, where the 
adiabatic condition 

3 uj eq e B , . 

- — V « 1 29 

holds (hence the unusual scaling with the charge, (eq e ) 8 ^ 3 ). Note that in this regime the 
process is independent of the MCP vacuum mass, m^, but this parameter still enters 
through the adiabatic condition ( |29p . For the value of k in the non-adiabatic regime we 
refer to Appendix A of |61j . 

The MCP pair production process damps the photon flux according to the usual 
decay law, so the photon survival probability after traveling physical distance L will be 
given by 

V(L) = exp (- J K(L)dLj , (30) 

where L is redshift-dependent. Using the redshift dependencies quoted in the discussion 
around Eq. (jTTl) . k redshiftdjjl as k(1 + z). This leads us to 

V{z)=exp{-Ky{z)) , (31) 

where y(z) is the comoving distance to the source. Note that this expression can 
be recovered from the ALP case, Eq. (|14p . in the A — > limit and substituting 
3P/(2L) -»■ k. 

As was noted in [TU1 [62] the above expression does not hold in what is probably the 
most interesting case, in which the MCPs arise from kinetic mixing. In that situation, 
photon to hidden photon oscillations also have to be taken into account and, most 

In the non-adiabatic regime this scaling is not valid. 
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Figure 6. Confidence levels (68% and 95%) on the fl m — k plane for the pure MCP 
model (left panel). Dark blue contours show constraints from SN data only, light 
blue from H{z) data, and solid line contours from joint SN+H(z). The right panel 
shows the corresponding constraints in the mini-charge vs MCP mass plane assuming 
B = InG. For comparison we show the previous SN dimming constraints [TU] (also for 
B = InG), CMB bounds [55] and the most sensitive purely laboratory experiments, 
light polarization [53], tests of the Coulomb's law [53] and accelerator cavities [55] . 



surprisingly, they tend to suppress the photon disappearance! In this scenario both 
photons and hidden photons get an effective mass from the magnetic-field-dressed MCP. 
However, the coupling of hidden photons to the MCP particles is much stronger than 
the corresponding coupling for photons, so the refractive effect (the effective mass) is 
always larger. The large mass of the hidden photon acts to suppress the mixing angle 
between photons and MCPs, in a similar manner a large plasma frequency (effective 
mass for the photon) suppresses the mixing between photons and ALPs in the previous 
section. The photon survival probability saturates at a value 



where x is the kinetic mixing between photons and hidden photons. The interested 
reader can find further details about these arguments in [62J. The value of \ is usually 
restricted to be smaller than 10 -3 since it has a radiative origin, so there is no foreseeable 
constraint on this scenario from cosmic opacity. Therefore, we must focus our attention 
on the pure MCP scenario. 

Fig. [6] (left) shows 1 and 2-a joint confidence levels in the k — fl m plane, again for 
SN data only (dark blue contours) and for the combined SN+if(2:) data set (solid line 
contours). In this case, SN data alone allow a zero cosmological constant in the presence 
of MCPs with a value k ~ 2.3 x 10 -4 as suggested in [TO]. However, the inclusion of 
H{z) data rules out this possibility and sets a strong bound 




(32) 




(33) 
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which, translated into MCP parameters, allows us to constrain the region in parameter 
space shown in Fig. [6] (right). 

5. Forecasts for future baryon acoustic oscillations and Supernovae surveys 

So far we have investigated constraints on cosmic opacity - and also their implications for 
models which violate photon number conservation - that are imposed from current data, 
namely from direct measurements of cosmic expansion H(z) using cosmic chronometers 
combined with Type la Supernova data (in particular the SCP Union 2008 compilation). 
However, new and more accurate data for H(z) (as well as dA_{z), the angular diameter 
distance) will be available through ongoing and future Baryon acoustic oscillations 
(BAO) surveys. In this section, we show forecasted constraints for cosmic opacity and 
the related models of sections [SHU that can be achieved by combining Supernova and 
H{z) data from future spectroscopic BAO surveys. We focus in particular on two BAO 
missions, namely the Baryon Oscillation Spectroscopic Survey (BOSS) and EUCLID. 
Finally, we also consider forecast constraints from proposed SN missions, in particular 
combining EUCLID and SNAP. 

BOSS [66] is part of the SDSS-III survey and is scheduled to operate over the 
period 2009-2014. Using the 2.5 m SDSS telescope, it will measure redshifts of 1.5 
million luminous galaxies in the range 0.1 < z < 0.7 (as well as Lya absorption towards 
160,000 high-redshift quasars at about z ~ 2.5), covering ~ 10, 000 deg 2 of high-latitude 
sky. The forecast precision for H(z) is 1.8%, 1.7% and 1.2% in redshifts bins centered 
at z — 0.35, 0.6 and 2.5 respectively. On the other hand, EUCLID - proposed to 
ESA's Cosmic Visions programme - aims for lunch around 2018. A combination of 
the earlier SPACE [67] and DUNE [68] missions, EUCLID would provide around 150 
million redshifts in the range z < 2, covering about 30,000 deg 2 . Fig. [7] shows forecast 
errors around the WMAP7 AC DM model for both BOSS and EUCLID. Also shown for 
comparison are the current H(z) "cosmic chronometers" data used above. 

We use the code developed by Seo & Eisenstein [69] to estimate the errors in radial 
distances achievable by using BAOs as a standard ruler. Fig. M shows our forecasted 
constraints on the parameter e of section [2], using the current type la SN data (Union 
2008) in combination with modelled BAO data with forecasted errors for both BOSS and 
EUCLID0. Note that although BOSS will achieve much smaller error bars than those of 
current H(z) data (cf Fig. [7]), it will span a much narrower redshift range, so it will in 
fact provide somewhat weaker constraints than the current H(z) "chronometers" data. 
To make a more direct comparison we have also shown the corresponding constraints 
obtained by restricting the current H(z) data in the narrower redshift range available 
to BOSS (thin solid line labeled "chronometer (low z)"). On the other hand, significant 
improvement of these constraints will be achieved by EUCLID. 

Similarly, in Figs. [9] & \IU\ we show forecasted constraints for the simple ALP 

+ We will consider forecasted constraints combining planned SN missions' data as well at the end of 
this section. 



Constraints on cosmic opacity and BSM physics from cosmological distances 



21 




o I ■■■■■■■■■■■■■■ ■ ■ ... I 
0.0 0.5 1.0 1.5 2.0 

z 

Figure 7. Forecasted error-bars on H{z) for BOSS (red) and EUCLID (green) 
compared to current H(z) data from differential ageing of galaxies (black). Also shown 
in the blue rectangle is the Ricss ct al. determination of Hq [24] . 

model of section |3] (A = 2/3) and for chameleons with A £ [2/3,4/3] (section |3]). 
For ALPs, EUCLID will provide significant improvement, notably by a factor of 2-3 
on the constraints on P/L in the weak mixing regime of Fig. El As the probability of 
mixing is inversely proportional to the square of the energy scale of the ALP coupling, 
this will result in a modest improvement of the bounds on M by factor of order unity. 
Finally, forecast constraints for MCPs (section H]) are displayed in Fig. [TU where, again, 
EUCLID will improve constraints on the parameter k by a factor of 2-3. As k ~ q^ 3 , 
this results in a modest, order unity improvement of the constraints on the charge of 
any exotic MCP. 

So far we have considered the effect that future BAO data will have on the 
constraints of sections EJHJ when combined with current SN data, and in particular 
the SCP Union 2008 compilation. We showed that EUCLID will lead to a significant 
improvement of these constraints (figures IHl fTTT) . while the narrower redshift range of 
BOSS renders it comparable to current H(z) measurements for constraining these 
models. We end this section by considering the effect that proposed SN surveys data will 
have on these constraints, in particular forecasted constraints for SNAP (or dark energy 
task force stage IV SNe mission) [70] combined with EUCLID. Fig. [12] summarises 
these constraints for all models considered above. Our forecasted constraints, shown 
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Figure 8. Forecasted constraints on the opacity parameter e of section[2]for BOSS and 
EUCLID, combined with current SN data. Also shown are current constraints from 
H(z) "chronometer" data (again joint with SN). The small shift towards negative e is 
due to the fact that actual SN data are being used. 

in orange scale, appear on top of the corresponding joint SN+H(z) constraints from 
current data, see sections 3 and 4. SN data from these proposed surveys will lead to 
notable improvement of these constraints, for example, by nearly an order of magnitude 
in the parameter e described above. As the figure shows, for the models considered 
above, this will correspond to an improvement of a factor of up to few on the strength 
of the coupling of ALPs to photons, and on the charge of MCPs. 

6. Conclusions 

If new particles from physics beyond the standard model couple to photons then the 
propagation of light may be altered. In this paper we have focused on two scenarios for 
exotic particles which can significantly modify the propagation of photons as they pass 
through magnetic fields. Measurements of cosmic opacity are a strong tool to constrain 
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Figure 9. The joint constraints of Fig. |3] for the simple ALP model of section |3] with 
A = 2/3, shown together with the corresponding forecast constraints from BOSS and 
EUCLID combined with current SN data. 




Figure 10. The joint constraints of Fig. [5] for chameleons with 2/3 < A < 4/3), 
together with the corresponding forecast constraints from BOSS and EUCLID joint 
with current SN data. 



such scenarios, as interactions between photons and exotic particles in the magnetic 
fields of the intergalactic medium leads to a new source of cosmic opacity. Uniform 
deviations from cosmic transparency (i.e. opacity) can be constrained through their 
effects on distance duality, by parameterizing possible deviations from the Etherington 
relation. The Etherington relation implies that, in a cosmology based on a metric theory 
of gravity, distance measures are unique: the luminosity distance is (1 + z) 2 times the 
angular diameter distance. Both luminosity distance and angular diameter distance 
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Figure 11. Constraints of Fig. |6]for MCPs, compared to forecasted constraints from 
BOSS and EUCLID, combined with current SN data. 
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Figure 12. Forecast constraints from joint SNAP+EUCLID (orange scale), shown 
together with the corresponding constraints from current data, namely SN (Union08) 
joint with chronometer H(z) (blue scale). Dark and light contours correspond to 1- and 
2-er respectively. From top left to bottom right: constraints on the opacity parameter 
e (section [2J , parameter P/L for the simple ALP model of section G2 parameters A & 
P/L for chameleons (section [3|) , and parameter k for MCPs (section 0|) . 
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depend on the Hubble parameter H(z), but this relation is valid in any cosmological 
background where photons travel on null geodesies and where, crucially, photon number 
is conserved. We have restricted our attention on violations of the Etherington relation 
arising from the violation of photon conservation. 

We have combined direct measurements of cosmic expansion (from the latest 
determinations of the Hubble parameter) at redshifts < z < 2 and recent SN 
data yielding the luminosity distance. SN-inferred luminosity distances are affected by 
violation of photon conservation, but the H(z) measurements we use are not. Assuming 
an underlying flat ACDM model, we have placed tight limits on possible deviations 
from photon-conservation. Photon-conservation can be violated by simple astrophysical 
effects which give uniform attenuation such as gray dust. We have reported updated 
constraints on this effect. 

More exotic sources of photon-conservation violation involve a coupling of photons 
to particles beyond the standard model of particle physics. We have focused on axion-like 
particles, new scalar or pseudo scalar fields which couple to the kinetic terms of photons, 
and mini-charged particles which are hidden sector particles with a tiny electric charge. 
Photons passing through intergalactic magnetic fields may be lost by pair production of 
light mini-charged particles. If the mixing between axion-like particles and photons is 
significant, then interactions in the intergalactic magnetic fields will also lead to a loss 
of photons due to conversion into ALPs. However if the coupling between photons and 
ALPs is sufficiently strong, one-third of any initial flux will be converted into ALPs, 
and two-thirds into photons, resulting in a redshift-independent dimming of supernovae 
which we cannot constrain or exclude with cosmic opacity bounds. 

The improved measurement of the cosmic opacity found here leads to improved 
bounds on these exotic physics scenarios which are summarised in Fig. [131 Future 
measurements of baryon acoustic oscillations, and an increase in the number of 
observations of high redshift supernovae will lead to further improvements in the 
constraints on physics beyond the standard model. 
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Figure 13. Prospects of the future bounds on ALPs (left) and MCPs (right). ALPs: 
The dark region corresponds to the current bounds while gray and light grey regions 
cover the possible future bounds using EUCLID and EUCLID+SNAP data. Taking 
into account the redshift dependence of the background improves these bounds in 
the region covered by the lines. MCPs: Current bounds (solid) and prospects from 
EUCLID (dashed) and EUCLID+SNAP (dotted). 



Appendix 

Here we summarise the relevant equations describing how the luminosity distance 
measures are affected by the presence of particles beyond the standard model that 
couple to photons: axion-like particles (ALP), chameleons and mini-charged particles 
(MCP). 

The presence of ALPs will have an impact on observations of supernovae if they 
are observed through magnetic fields. The strength and spatial structure of the inter- 
galactic magnetic fields is highly unconstrained and so they are left as parameters. The 
luminosity distance to SNe is modified by the redshift-dependent probability of photon 
survival V{z): 

d L (z)^d L {z)/^/V{z). (34) 

There are two physical limits in which V(z) can be analytically calculated, namely 
the incoherent and coherent regimes. In the incoherent regime one has: 

V(z) = A + (1 - A) exp {-\^P) , (35) 

where L is the size of magnetic domains, P the probability of photon to ALP conversion, 
and y(z) the comoving distance to the source. In the coherent regime, one can take 
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into account additional effects due to evolution of the background magnetic fields with 
redshift. In this case the probability of photon survival reads: 

V{z) = A + (1 - A) exp (-jfj^^) , (36) 

where H and Q m are the Hubble constant and matter density parameter (at the present 
epoch) respectively, and H(z) is the Hubble parameter at redshift z. For the simplest 
ALP models, A = 2/3 in the above equations, while for chameleons these equations 
hold with A ^ 2/3. 

The existence of low-mass MCPs can also have a tremendous impact on photon 
propagation over cosmological distances. Photons from a given source can for instance 
pair produce MCPs, even without the need for a (CMB) photon, in the presence of an 
inter-galactic magnetic field. The MCP pair production process damps the photon flux, 
so, again, the luminosity distance to SNe is modified by a redshift-dependent probability 
of photon survival, V(z). In this case one has: 

P(z)=eKp(-Ky{z)) , (37) 

where k is the photon to MCP transition rate and is given by 

Here, B denotes the magnetic field strength, u is the photon energy, q e is the MCP 
electric charge in units of the electron charge e, and / is an order unity factor which 
depends on the nature of the MCP and the photon polarization with respect to the 
magnetic field. T denotes the usual T-function. 



/ • (38) 
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